Dynamic melting of confined vortex matter 
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We study dynamic melting of confined vortex matter moving in disordered, mesoscopic channels 
by mode-locking experiments. The dynamic melting transition, characterized by a collapse of the 
mode-locking effect, strongly depends on the frequency, i.e. on the average velocity of the vortices. 
The associated dynamic ordering velocity diverges upon approaching the equilibrium melting line 
Tm,e{B) as Vc ~ {Tm,e — T)~^ . The data provide the first direct evidence for velocity dependent 
melting and show that the phenomenon also takes place in a system under disordered confinement. 
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An intriguing aspect of periodic media driven through 
a (random) pinning potential, e.g. vortex lattice (VL's) 
in superconductors [1], charge density waves (CDW's) 
[2] , or conventional sliding solids [3] , is the possibility of 
a dynamic ordering (DO) transition [4,5]: while at small 
velocity pinning may disrupt the lattice, causing liquid- 
like flow [6] around pinned islands, at large velocity the 
influence of pinning diminishes and the elastic interac- 
tions dominate, favoring a crystalline structure. 

The DO phenomenon becomes particularly interesting 
close to the thermodynamic melting point of the medium, 
since in addition to fluctuations due to pinning also ther- 
mal fluctuations become relevant [5,7]. Their combined 
effect was first studied by Koshelev and Vinokur [5] for 
2D VL's. They introduced the concept of a shaking tem- 
perature Tsh, characterizing the fluctuations in the mov- 
ing frame due to the quenched disorder. This shaking 
diminishes with increasing velocity as ksTsu = ^p,v/v 
(Tp^^ reflects the pinning strength and viscous damping) . 
Fluid-like, incoherent motion sets in when the effective 
'temperature' T+T^h equals the equilibrium melting tem- 
perature Tm,e of the pure system. Hence, for T < Tm,e 
there exists a recrystallization velocity 



^p,v/[kB{T„ 



(1) 



which diverges when temperature approaches Tm,e- Indi- 
cations for this scenario were found in experiments on su- 
perconducting fllms [8] by identifying an inflection point 
in the current-voltage {IV) curves with DO. This identi- 
fication seemed justified from simulations [9] but more di- 
rect experimental evidence for the phenomenology in [5] 
is still lacking. Simulations of a density gradient driven 
system [10] in fact suggest that the inflection point can 
also be due to changes in large-scale flow morphology 
while recent experiments on crystals [11] have shown that 
for a 3D VL the inflection point is due to macroscopic co- 
existence of two phases and not due to microscopic DO. 

In this Letter we study DO of confined vortex matter 
flowing in disordered mesoscopic channels using a com- 
pletely different dynamic probe. In [12] we showed how a 
mode-locking (ML) technique can be used to explore the 



flow conflguration in the channels in presence of strong 
disorder from the vortex arrays in the channel edges. The 
ML phenomenon occurs due to coupling between, on the 
one hand, lattice modes of frequency fint = Qv/a, with 
q an integer and a the lattice periodicity, which may ex- 
ist in an array that moves coherently with velocity v, 
and, on the other hand, a superimposed rf-drive of fre- 
quency / at an integer fraction 1/p of fint [13]. This cou- 
pling produces plateaus in the dc-transport (IV) curves 
or sharp peaks in the differential conductance (dl/dV) 
when V — {p/q)fa. However, these peaks are reduced 
and can eventually vanish due to incoherent velocity fluc- 
tuations in the moving array, arising from both thermal 
and quenched disorder [14,15]. Thus, the collapse of the 
ML peak marks the transition from coherent to incoher- 
ent flow. Using this criterion we provide for the flrst 
time conclusive evidence for a velocity dependent melt- 
ing transition at Vc [16] and probe the divergence of Vc 
upon approaching the static phase boundary. 

The experiments are performed on a superconducting 
double layer of weak pinning amorphous (a-)Nbi_a;Gea; 
{x ~ 0.3, 550 nm thickness, Tc = 2.68 K, normal re- 
sistivity pn ~ 2/if2m) and strong pinning NbN (50 nm 
thickness) on top, containing Nch{— 200) straight chan- 
nels of width w = 230 nm etched to a depth of 300 nm 
[17]. Vortices in the channels are conflned between chan- 
nel edges (CE's) consisting of pinned, disordered VL's 
which impose both a periodic and random potential on 
the conflned vortices via their mutual shear interaction, 
characterized by a shear modulus cge- This potential de- 
termines both the threshold force density Fp ~ cee/w 
and the dynamics of vortices in the channel [12,18,19]. 
We measured dc and dc-rf transport versus magnetic fleld 
and temperature using a four probe conflguration. The 
sample was immersed in superfluid ^He. The frequency 
of the applied rf current was in the range 1 — 200 MHz, 
its amplitude Irf could be as large as 2 mA [20]. 

Due to the confinement, the ML condition for chan- 
nel fiow attains a particularly useful form: although the 
array in the channel may be frustrated, meaning that 
the average longitudinal periodicity a ag and the row 



1 



spacing b ^ bo (ao = 26o/\/3 ~ 1.075^/ ^q/B arc the 
equilibrium values), the voltage at which the main inter- 
ference {p = q = 1) occurs is simply given by [12] 

= $o/niVe^, (2) 

with n the number of moving rows in each channel. 




B(T) 

FIG. 1. Shear force density (determined using a veloc- 
ity criterion v/ao ~ 1 MHz) versus field at a temperature 
T = 1.94 K. Vertical lines indicate the field of the transition 
from n = 6 — > 7 and n = 7 — *■ 8, dashed vertical line marks 
the field where the n = 8 — > 9 transition is expected. Open 
symbols (right axis): field dependence of the zero bias resis- 
tance Rq. Inset: dl /dV versus normalized voltage in presence 
of 85 MHz rf-currents for B = 0.9 - 1.16 T (steps of 20 mT). 
The data illustrate the n = 7 — »■ 8 transition in the ML signal. 

We first discuss the typical behavior in the solid phase 
where matching cfFcxits between the confined array and 
the channel width dominate the behavior. In Fig. 1 the 
thick line shows the dc-depinning force density Fp = JpB 
[Irf = 0). For fields _B < 1.1 T two oscillations in Fp are 
seen. The inset shows dl/dV curves in presence of an 
85 MHz rf-current (7^/ = 0.53 mA) versus voltage for 
fields 0.9 T< B < 1.16 T. The data for 0.9 T exhibits a 
peak in dl/dV corresponding to ML of 7 vortex rows in 
each channel. For larger fields the amplitude of this peak 
decays while another peak, corresponding to 8 rows, ap- 
pears. The field region where these peaks coexist (0.92 
T ^ B < 1 T) evidently corresponds to the situation of 
maximum mismatch between the natural width nbo of 
n = 7, 8 vortex rows and the effective channel width (es- 
timated as Weff = 315 nm). In the same field regime Fp 
exhibits a maximum. As shown in [12], this maximum at 
mismatch is caused by jamming of the flow at locations 
in the channel where the number of rows switches from n 
to n± 1. The motion there is partially blocked by disloca- 
tions with Burgers vectors that are almost perpendicular 



to the fiow direction. The structural disorder of the fixed 
VL in the CE's is responsible for this phenomenon. 

We now turn to the behavior in larger magnetic fields, 
S > 1.1 T. At B ~ 1.15 T, Fp shows a rapid decrease 
after a small upturn. This drop coincides with the on- 
set of a measurable zero-bias resistance Rq (displayed on 
the right axis). These two distinct features reflect the 
loss of shear rigidity of the vortex configuration inside 
the channel, indicating a transition to a confined vortex 
liquid [17]. Note that these features occur well below the 
magnetic field at which the transition from n = 8 to 9 
vortex rows is expected from the condition 8.56o = Weff, 
namely Bg^g = 1.27 T. In addition, neither Fp nor Rq 
show any particular features at i?8,9) which shows that 
the liquid is insensitive to (mis)matching eflFects. 




FIG. 2. Differential conductance versus reduced voltage at 
T = 1.94 K in fields ranging from 1.14 T (topcurve) to 1.3 T 
(20 mT steps) for (a) rf-current of 6 MHz (7^/ = 0.09 mA) 
and (b) 140 MHz and Irf = 0.94 mA. 

The ML experiments provide important new informa- 
tion regarding the coherence and the shear rigidity of the 
moving array. Figure 2(a) shows dl/dV data in the field 
range 1.14 — 1.3 T plotted versus the normalized voltage 
for a 6 MHz rf-current. The ML peak in the upper curve 
(7? = 1.14 T) corresponds to the coherent motion of 8 
rows. For larger field the peak amplitude drops rapidly 
and vanishes at B ~ 1.2 T. Above 1.2 T the dl/dV ver- 
sus V curves remain featureless. The vanishing of the 
ML signal marks the transition from coherent motion of 
a vortex solid to incoherent 'liquid' motion. A similar 
interpretation of this phenomenon was suggested in [14] 
with respect to the VL melting transition in YBCO. 

In Fig. 2(b) a dl/dV data set for the same field range 
is shown, but now for a 140 MHz rf-current. The average 
vortex velocity at mode-locking, v — fa is therefore over 
20 times larger than in Fig. 2(a). Again a clear ML peak 
is observed at B = 1.14 T, but in contrast to the 6 MHz 
data in (a) the field range in which ML takes place is 
considerably larger, extending up to B ~ 1.3 T. The na- 
ture of the moving medium thus depends strongly on its 
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velocity. For example: at B = 1.24 T one observes 'liq- 
uid' motion for / = 6 MHz but coherent motion at 140 
Mhz. Interestingly, in the early work of Fiory [13] it was 
already mentioned that the vanishing of the ML signal 
shifts to larger fields when measured at larger frequency. 
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FIG. 3. (a) ML step width A7, as defined in the inset to 
(b), versus field for / = 6 MHz and / = 60 MHz and several 
values of Irf- Arrows indicate the dynamic field Bc{f). (b) 
Data points: fc versus B. Drawn line is a fit according to 
fc ~ {Bm,e — B)^^ . The dashed lined marks Bm,e = 1-36 T. 

For a proper determination of the collapse of the ML 
signal we should take into account its dependence on Irf 
(see e.g. [21]). Therefore the ML step width A/ (defined 
in the inset to Fig. 3(b)) was determined as fmiction of 
field for various rf-amplitudes and frequencies. The re- 
sult is shown in Fig. 3(a) for 6 MHz and 60 Mhz. As 
observed, A/ vanishes linearly with _B at a dynamic melt- 
ing field Bc{f) which is essentially independent of 7^/. 

Next we study the frequency dependence of B,,- Rather 
than plotting versus /, we present the results by plot- 
ting the frequency at which ML vanishes, denoted by /c, 
versus i?, as done in Fig. 3(b). At low fields fc slowly 
increases with field, but it starts to diverge at larger field. 
Such divergence is expected when we identify Vc = fed as 
the ordering velocity proposed in [5] . In the present case, 
the plastic (or 'liquid') flow for v < Vc is controlled by 
the disordered VL's in the CE's, as we concluded from 
simulations at small drive (as in [12,19]). This plastic 
flow involves frequent slip events and transverse vortex 
jumps between rows. Such behavior has also been ob- 
served in magnetic bubble arrays moving along a rough 
wall [22]. For v > Vc these eflFects disappear and co- 
herent motion sets in. For B > 1.32 T we could no 
longer resolve the ML effect below 200 MHz. In this fleld 



regime thermal fluctuations alone are sufficient to induce 
incoherent motion, regardless of the dynamic influence 
of disorder. The fc{B) data can approximately be fitted 
to fc ^ {Brn.e ^ B)~'^'^ with 1/^ = 2 and an equilibrium 
melting field Bm,,e = 1-36 T. An exponent vb = 1, as 
expected from Eq.(l) and the behavior of the melting 
line, yields a rather poor fit. This discrepancy with the 
data may originate from the fact that with changing field 
not only the 'distance' Bm,e — B to the phase boundary 
varies, but also the commensurability, i.e. the dislocation 
density and effective disorder in the channel. 
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FIG. 4. (a) Data points: dynamic ordering frequency at 
B = 1.16 T versus temperature. Drawn line: fit according 

to /, ~ (T,-n.e - T)-^ with Tm,e = 2.011 K (indicated by the 
dashed line), (b) Dc-resistance V/I {Irf = 0) versus T for 
several dc-currents I. 

To avoid the possible influence of mismatch and pro- 
vide a more direct test of Eq.(l), we determined fc as a 
function of temperature for i? = 1.16 T (near matching) 
from an analysis of the ML amplitude A/ at several rf- 
currents similar to that in Fig. 3(a) (details are given in 
[23]). The result is shown in Fig. 4(a). Again a clear di- 
vergence of fc is observed. In this case the data can be fit- 
ted quite well to Eq.(l) written as fc = fo{^ — T/Tm,e)~^ 
with /o = 0.174 MHz and T^^e = 2.011 K. Interestingly, 
a comparison with the dc-resistance at various currents 
(Fig. 4(b)) shows that the equilibrium melting tempera- 
ture Tm,e coincides with the temperature where the R{T) 
curves merge, i.e. the IV curves become fully linear. 
This confirms the usual assumption made in dc-transport 
studies on VL melting [24]. In measurements as func- 
tion of field we observed the merging of Ri<.i^a{B) and 
Ri=iooha{B) at S ~ 1.34 T, in good agreement with 
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Bm,e determined above. We further note that none of 
the dc-IV curves show an inflection point, implying that 
such feature is not required for microscopic DO [11,25]. 

It is remarkable that our system, in which 'edge disor- 
der' [26] is the prime source of pinning, shows DO typical 
for a VL with bulk disorder. Let us therefore discuss the 
nature of the DO in more detail. Even for a 2D VL with 
bulk pinning the DO was under intensive debate [27-30]. 
At present, the most likely scenario seems that, on in- 
creasing V, after a crossover from fully plastic to partially 
layered, smectic flow, finally a transition to a Moving 
Transverse Solid occurs [28]. At this transverse freezing 
transition (TFT), inter-chain excursions (so called per- 
meation modes [29]) are suppressed, but free dislocations 
with Burgers vector parallel to iJ [30] and only longitu- 
dinal short-range order (LSRO) remain. Particularly, it 
is the TFT which is described by T^h [28]. Turning to 
the channels, we then suggest that our DO boundary 
Vc{T,B) reflects the TFT in which permeation modes 
due to roughness of the CE arrays are suppressed. Pre- 
liminary simulations support this view [23] . Additionally, 
above Vc we observe only incom,plete ML, suggesting in- 
deed LSRO [29] and residual slip between chains. 

The shaking effect in the channels is estimated as fol- 
lows. The characteristic frequency /o — p.v I [kBTm.eo) 
obtained from Eq.(l), can be derived from Ref. [5] as: 
/o = 'i / iiT^uP J / {'^^a^ dksTrn.e) with pj- the flux flow 
resistivity, d the film thickness and 7„ the pinning en- 
ergy sqiiarcd times the 2D pinning range. For the chan- 
nels, the short wavelength oq) disorder component 
due to vortex displacements d in the CE acts in a range 
~ ao/2 from the CE's and has a strength ~ AgCge with 
A, ~ (i/(|d|2)/ao)/(7rV3) [19,23]. We thus assume that 
'shaking' of the first vortex layer near each CE domi- 
nates the TFT. Hence, 7„ ~ (AsC66ao&orf)^(ao/2)^. Us- 
ing the melting criterion ATikBTm.e — c^^a^d one ob- 
tains /o ~ 2QAlpfkBTm.el^ld. Taking pf ~ /9„/2 and 
d = 300 nm yields /o = • 500 MHz. This is in rea- 
sonable agreement with the measured value /o = 0.174 
MHz when wc assume Ag ~ 0.02, i.e. rms relative dis- 
placements in the CE of {\d\^)/ oq) ~ 0.1. 

In conclusion, dynamic melting of vortex matter driven 
through disordered channels was studied by mode-locking 
experiments. The melting line strongly depends on the 
ML frequency, i.e. the average velocity. The associated 
ordering velocity diverges upon approaching the equilib- 
rium melting line, yielding a dynamic phase diagram with 
coherent, plastic and fiuid fiow as predicted theoretically 
[5,7] . The ML technique presents a powerful tool to study 
phase transitions in driven periodic media and we hope 
our results will stimulate similar investigations in related 
fields like CDW dynamics and solid friction. 
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